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Abstract — The availabihty of positional information is of great 
importance in many commercial, public safety, and military 
applications. The coming years will see the emergence of location- 
aware networks with sub-meter accuracy, relying on accurate 
range measurements provided by wide bandwidth transmissions. 
In this two-part paper, we determine the fundamental limits of 
localization accuracy of wideband wireless networks in harsh 
multipath environments. We first develop a general framework 
to characterize the localization accuracy of a given node here and 
then extend our analysis to cooperative location-aware networks 
in Part II. 

In this paper, we characterize localization accuracy in terms 
of a performance measure called the squared position error 
bound (SPEB), and introduce the notion of equivalent Fisher 
information to derive the SPEB in a succinct expression. This 
methodology provides insights into the essence of the localization 
problem by unifying localization information from individual 
anchors and information from a priori knowledge of the agent's 
position in a canonical form. Our analysis begins with the re- 
ceived waveforms themselves rather than utilizing only the signal 
metrics extracted from these waveforms, such as time-of-arrival 
and received signal strength. Hence, our framework exploits all 
the information inherent in the received waveforms, and the 
resulting SPEB serves as a fundamental limit of localization 
accuracy. 

Index Terms — Cramer-Rao bound (CRB), equivalent Fisher 
information (EFI), information inequality, localization, ranging 
information (RI), squared position error bound (SPEB). 



I. Introduction 

Location-awareness plays a crucial role in many wireless 
network applications, such as localization services in next 
generation cellular networks IT], search-and-rescue operations 
lUl, El, logistics f4l, and blue force tracking in battlefields fS). 
The Global Positioning System (GPS) is the most important 
technology to provide location-awareness around the globe 
through a constellation of at least 24 satellites f^"], [7|. How- 
ever, the effectiveness of GPS is limited in harsh environments, 
such as in buildings, in urban canyons, under tree canopies, 
and in caves |8|, |9|, due to the inability of GPS signals to 
penetrate most obstacles. Hence, new localization techniques 
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Fig. 1. Location-aware networks: the anchors (A, B, C, and D) communicate 
with the agents (1 and 2), and each edge denotes a connection link between 
anchor and agent. 



are required to meet the increasing need for accurate localiza- 
tion in such harsh environments |8|, |9|. 

Wideband wireless networks are capable of providing accu- 
rate localization in GPS-denied environments ll8l- lfT2l . Wide 
bandwidth or ultra-wide bandwidth (UWB) signals are par- 
ticularly well-suited for localization, since they can provide 
accurate and reliable range (distance) measurements due to 
their fine delay resolution and robustness in harsh environ- 
ments |fT3l - ||20l . For more information about UWB, we refer 
the reader to llITl-lfSSl. 

Location-aware networks generally consist of two kinds of 
nodes: anchors and agents. Anchors have known positions 
(for example, through GPS or system design), while agents 
have unknown positions and attempt to determine their po- 
sitions (see Fig. [U. Each node is equipped with a wideband 
transceiver, and localization is accomplished through the use 
of radio communications between agents and their neighboring 
anchors. Localizing an agent requires a number of signals 
transmitted from the anchors, and the relative position of the 
agent can be inferred from these received waveforms using 
a variety of signal metrics. Commonly used signal metrics 
include time-of-arrival (TOA) 18], 0, |[l7|-E0l, Il27l-ll30l. 
time-difference-of-arrival (TDOA) |I3T1, [l32l, angle-of-arrival 
(AOA) IS, 133], and received signal strength (RSS) l|9). Il34l. 
[35|. 

Time-based metrics, TOA and TDOA, are obtained by 
measuring the signal propagation time between nodes. In ideal 
scenarios, the estimated distance equals the product of the 
known propagation speed and the measured signal propagation 
time. The TOA metric gives possible positions of an agent on 
a circle with the anchor at the center, and it can be obtained by 
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either the one-way time-of-flight of a signal in a synchronized 
network HU, HH, lH, iH, or the round-trip time-of-flight 
in a non-synchronized network ||26l . |[36|. Ahernatively, the 
TDOA metric provides possible positions of an agent on the 
hyperbola determined by the difference in the TOAs from two 
anchors located at the foci. Note that TDOA techniques require 
synchronization among anchors but not necessarily with the 
agent. 

Another signal metric is AOA, the angle at which a signal 
arrives at the agent. The AOA metric can be obtained using an 
array of antennas, based on the signals' TOA at each antennaQ 
The use of AOA for localization has been investigated, and 
many hybrid systems have been proposed, including hybrid 
TOA/AOA systems f30l, f^Tl, and hybrid TDOA/AOA sys- 
tems |42|. However, some of these studies employ narrowband 
signal models, which are not applicable for wideband antenna 
arrays. Others are restricted to far-field scenarios or use far- 
field assumptions. 

RSS is also a useful metric that uses the strength of the 
received signal to estimate the propagation distance between 
nodes Q, If34l . ||36l . This metric can be measured during the 
data communications using low-complexity circuits. Although 
widely implemented, RSS has limited accuracy due to the 
difficulty in precisely modeling the relationship between the 
RSS and the propagation distance E), jg). 

Note that the signal metrics extracted from the received 
waveforms may discard relevant information for localization. 
Moreover, models for the signal metrics depend heavily on 
the specific measurement processes0 Therefore, in deriving 
the fundamental limits of localization accuracy, it is necessary 
to utilize the received waveforms rather than the signal metrics 
extracted from the waveforms f28l, f29l, f46l, pT). 

Since the received waveforms are affected by random 
phenomena such as noise, fading, shadowing, multipath, and 
non-line-of-sight (NLOS) propagations |48|, |49|, the agents' 
position estimates are subject to uncertainty. The Cramer- 
Rao bound (CRB) sets a lower bound on the variance of 
estimates for deterministic parameters f5Q\, fST], and it has 
been used as a performance measure for localization accuracy 
||52l . However, relatively few studies have investigated the 
effect of multipath and NLOS propagations on localization 
accuracy. Multipath refers to a propagation phenomenon in 
which a transmitted signal reaches the receive antenna via 
multiple paths. The superposition of these arriving paths 
results in fading and interference. NLOS propagations, created 
by physical obstructions in the direct path, produce a positive 
bias in the propagation time and decrease the strength of the 
received signal, which can severely degrade the localization 
accuracy. Several types of methods have been proposed to deal 
with NLOS propagations: 1) treat NLOS biases as additive 

'The AOA metric can be obtained in two ways, directly through measure- 
ment by a directional antenna, or indirectly through TOA measurements using 
an antenna array ( 37 |-| 40 |. Wideband directional antennas that satisfy size and 
cost requirements are difficult to implement, since they are required to perform 
across a large bandwidth |36|. As such, antenna arrays are more commonly 
used when angle measurement for wide bandwidth signals is necessary. 

-For instance, the error of the TOA metric is commonly modeled as an 
additive Gaussian random variable [8J, [30] , 143J. This model contradicts the 
studies in [T^-(2^, (44), (45|, and the experimental results in (8), fT6l . 



noise injected in the true propagation distances (SJ, Il53l . ll54l FI 
2) identify and weigh the importance of NLOS signals for 
localization [55J- [60J ; or 3) consider NLOS biases as param- 
eters to be estimated ||27l-||30l, ||46l. El, ED, ISlil- The 
authors in fSl, f^l, ||28l, ||29l showed that NLOS signals do not 
improve localization accuracy unless a priori knowledge of the 
NLOS biases is available, but their results were restricted to 
specific models or approximations. Moreover, detailed effects 
of multipath propagations on localization accuracy remains 
under-explored. 

In this paper, we develop a general framework to determine 
the localization accuracy of wireless networks^ Our analysis 
begins with the received waveforms themselves rather than 
utilizing only signal metrics extracted from the waveforms, 
such as TOA, TDOA, AOA, and RSS. The main contributions 
of this paper are as follows: 

> We derive the fundamental limits of localization ac- 
curacy for wideband wireless networks, in terms of a 
performance measure called the squared position error 
bound (SPEB), in the presence of multipath and NLOS 
propagation. 

• We propose the notion of equivalent Fisher information 
(EFT) to derive the agent's localization information. This 
approach unifies such information from different anchors 
in a canonical form as a weighed sum of the direction ma- 
trix associated with individual anchors with the weights 
characterizing the information intensity. 

« We quantify the contribution of the a priori knowledge 
of the channel parameters and agent's position to the 
agent's localization information, and show that NLOS 
components can be beneficial when a priori channel 
knowledge is available. 

• We derive the performance limits for localization systems 
employing wideband antenna arrays. The AOA metric 
obtained from antenna arrays are shown not to further 
improve the localization accuracy beyond that provided 
by TOA metric alone. 

> We quantify the effect of clock asynchronism between 
anchors and agents on localization accuracy for networks 
where nodes employ a single antenna or an array of 
antennas. 

The rest of the paper is organized as follows. Section HI] 
presents the system model, the notion of the SPEB, and the 
Fisher information matrix (FIM) for the SPEB. In Section HIH 
we introduce the notion of EFI and show how it can help 
the derivation of the SPEB. In Section IIVI we investigate 
the performance of localization systems employing wideband 
antenna arrays. Section |V] investigates the effect of clock 
asynchronism between anchors and agents. Discussions are 
provided in Section |VI] Finally, numerical illustrations are 
given in Section IVIII and conclusions are drawn in the last 
section. 

^In practice, however, a NLOS induced range bias can be as much as a few 
kilometers depending on the propagation environment |48|, |55|, and small 
perturbation may not compensate for NLOS induced en'or 

''in Part II Ii63il , we extend our analysis to cooperative location-aware 
networks. 
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Notations: The notation Ezj-} is the expectation operator 
with respect to the random vector z; A >^ B denotes that the 
matrix A — B is positive semi-definite; tr{-} is the trace of a 
square matrix; [ • denotes the upper left n x n submatrix 
of its argument; [ • ]„ .,„ is the element at the nth row and mth 
column of its argument; || • || is the Euclidean norm of its 
argument; and the superscripts [ • ]^ represents the transpose 
of its argument. We denote by /(x) the probability density 
function (PDF) /x (x) of the random vector X unless specified 
otherwise, and we also use in the paper the following function 
for the FIM: 



F,(w;x,y) 



|-ln/(w) 
ox 



|-ln/(w) 



where w can be either a vector or a symbol^ 

II. System Model 

In this section, we describe the wideband channel model 
Oil, 123, 1261, IISl, formulate the problem, and briefly 
review the information inequality and Fisher information. We 
also introduce the squared position error bound, which is a 
fundamental limit of localization accuracy. 

A. Signal Model 

Consider a wireless network consisting of anchors and 
multiple agents. Anchors have perfect knowledge of their 
positions, and each agent attempts to estimate its position 
based on the received waveforms from neighboring anchors 
(see Fig. [TJ0 Wideband signals traveling from anchors to 
agents are subject to multipath propagation. 

Let p G denote the position of the agentj^ which is 
to be estimated. The set of anchors is denoted by A/b = 
{1, 2, . . . ,iVb} = A/l U A/nl, where TVl denotes the set of 
anchors that provide LOS signals to the agent and A/nl denotes 
the set of remaining anchors that provide NLOS signals to the 
agent. The position of anchor k is known and denoted by 
Pfc G (fc G A/i,). Let (pk denote the angle from anchor k to 
the agent, i.e.. 



tan 



-1 y-yk 

X ~ Xk 



where p = [x and p^ = [xk ykY ■ 

The received waveform at the agent from anchor k can be 
written as 



rk( 



,(t) = £4"s(t-rf ) tG[0,Tob), (1) 

1=1 

where s{t) is a known wideband waveform whose Fourier 
transform is denoted by S{f), af^ and t^/' are the amplitude 
and delay, respectively, of the llh path, Lk is the number 



'For example, w is replaced by symbol r\0 in the case that /(■) is a 
conditional PDF of r given 0. 

^Agents estimate their positions independently, and hence without loss of 
generality, our analysis focuses on one agent. 

'We first focus on two dimensional cases and then extend the results to 
three-dimensional cases where p G R'^. 



of multipath components (MFCs), Zk{t) represents the obser- 
vation noise modeled as additive white Gaussian processes 
with two-side power spectral density Nq/2, and [0,To[,) is 
the observation interval. The relationship between the agent's 
position and the delays of the propagation paths is 

1 



.(0 



Pfel 



u(0 



(2) 



where c is the propagation speed of the signal, and b'l^ > 
is a range bias. The range bias b'j^^ = for LOS propagation, 
whereas 6^''' > for NLOS propagation!! 



B. Error Bounds on Position Estimation 

Our analysis is based on the received waveforms given by 
([T]i, and hence the parameter vector 9 includes the agent's 
position and the nuisance multipath parameters l9|, l62l . i.e.. 



6 



where Hk is the vector of the multipath parameters associated 
with rk{t), given by 



4^^ 



1,(2) „,(2) 



-'k 



(1) ^(2) ^,(2) 



a 



^k 



k 



fee A/L, 

n T 

"fe J ' 
k G A/'nl ■ 



Note that b^j^' — for k G A/L and is excluded from Hk- 

We introduce r as the vector representation of all the 
received waveforms rk{t), given by 

r^[rj rl ••• , 

where is obtained from the Karhunen-Loeve expansion of 
fkit) l50l . BTl . Let 6 denote an estimate of the parameter 
vector based on observation r. The mean squared error 
(MSB) matrix of 6 satisfies the information inequality l50l . 

El, m 

Er,e{(0-0)(0-0)''} ^ Je'- (3) 

where 3g is the Fisher information matrix (FIM) for the 
parameter vector 00 Let p be an estimate of the agent's 
position, and it follows from (O tha0 



2x2 ■ 



*LOS propagation does not introduce a range bias because there is an 
unblocked direct path. NLOS propagation introduces a positive range bias 
because such signals either reflect off objects or penetrate through obstacles. 
In this paper, received signals whose first path undergoes LOS propagation are 
referred to as LOS signals, otherwise these signals are referred to as NLOS 
signals. 

'when a subset of parameters is random, Jg is called the Bayesian 
information matrix. Inequality jsj also holds under some regularity conditions 
and provides lower bound on the MSE matrix of any unbiased estimates of 
the deterministic parameters and any estimates of the random parameters |50|, 
|65|. With a slight abuse of notation, g{-} will be used for deterministic, 
hybiid, Bayesian cases with the understanding that the expectation operation 
is not performed over the deterministic components of 0. 

'"Note that for three-dimensional localization, we need to consider a 3 X 3 
matrix [Jg^]^^^. 
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and hence 



E, 



,.{||p-pin>tr{[j,i],^J 



(4) 



Therefore, we define the right hand side of (|4]i as a measure 
to characterize the Umits of position accuracy as follows. 

Definition 1 ( Squared Position Error Bound): The squared 
position error bound (SPEB) is defined to be 



^(P)^tr{[J-],,,} 



C. Fisher Information Matrix 

In this section, we derive the FIM for both deterministic 
and random parameter estimation to evaluate the SPEB. 

1) FIM without A Priori Knowledge: The FIM for the 
deterministic parameter vector is given by 1501 



Je =Fr(r|0;0,0)^ 



(5) 



where f{v\9) is the likelihood ratio of the random vector r 
conditioned on 9. Since the received waveforms from different 
anchors are independent, the likelihood ratio can be written as 

m 



fir\e)= n /(rfel^) 

fceA/i 



(6) 



where 



f{rk\9) ccexpi — 



Tub 



1=1 



L;=i 



No Jo 

Substituting Q in we have the FIM Jg as 



dt 



1 



>e = 



TlAlTl + TnlAnlTjsjl 
AnlTnl 



Tnl-Anl 
Anl 



(7) 



where Al, Tl, Anl, and Tnl are given by dTTT ) and (l42b . In 
the above matrices, Al and Tl are related to the LOS signals, 
and Anl and Tnl are related to the NLOS signals. 

2) FIM with A Priori Knowledge: We now incorporate 
the a priori knowledge of the agent's position and channel 



The joint PDF of observation and parameters can be written 

as 

f{r,e) = f{r\e)f{e), 

where f{r\9) is given by and hence the FIM becomes 

(10) 



'p ' 



where Jw = Fr,e(r|6» ; 6», 6/) and Jp = Fe(6/;6/,6») are 
the FIMs from the observations and the a priori knowledge, 
respectively^ The FIM can be obtained by taking the 
expectation of Je in ^ over the random parameter vector 9, 
and Jp can be obtained by substituting ^ in ( fTOl ) as 



3p + E 



jfc 



' — *K,K 







(11) 



where Hp describes the FIM from the a priori knowledge of 
p, given by 

Hp = Fe(p;p,p) , 

and H^_^ = Fg{Kk\p ; Kk, Kk), H^^p = Fe(/tfc|p ; p, p), 
and Hp^ = Fe(/tfe|p ; p, /t^) characterize the joint a priori 
knowledge of p and Kk. 

D. Equivalent Fisher Information Matrix 

Determining the SPEB requires inverting the FIM Jg in 
(|7]l and ( fTOb . However, Je is a matrix of high dimensions, 
while only a small submatrix [Je^]2x2 interest. To 

circumvent direction matrix inversion and gain insights into 
the localization problem, we first introduce the notions of EFI 



Definition 2 (Equivalent Fisher Information Matrix): 
Given a parameter 9 — '[6\ and the FIM Je of the 

form 



A 



B 

C 



(12) 



parameters for localization. Since the multipath parameters Kh information matrix (EFIM) for 9i is given b 



are independent a priori, the PDF of 9 can be expressed aa"l 



where 9 e M^, 9i e M", A e M"^", B e 

and C e ^{N-n)x(N-n) ^jjjj „ < jj^g equivalent Fisher 

30 



/w=/(p) n /('^'^ip) 

fceA^b 



Je(6/i) = A-BC-^B^ 



(13) 



(8) 



where /(p) is the PDF of the agent's position, and /(ka;|p) 
is the joint PDF of the multipath parameter vector Kk con- 
ditioned on the agent's position. Based on the models of 
wideband channels ||36l, gOl, El and UWB channels (Hi, 
lITlJ . [24], |26l, p36l, we derive /(Kfe|p) in ^ in Appendix 
Island show that 



(9) 



where 



fe 



P/cl 



"when a subset of parameters are deterministic, they are eUminated from 



Note that the EFIM retains all the necessary information 
to derive the information inequality for the parameter vector 
9i, since [Je^]„xn = JsTH^OO and the MSE matrix of 
the estimates for 9i is bounded below by 3^^{9i)^ For 
two-dimensional localization {n = 2), we aim to reduce the 
dimension of the original FIM to the 2 x 2 EFIM. 

'^Note that Jg in IIP) requires averaging over the random parameters, and 
hence does not depend on any pailicular value of 0. In contrast, Jg in (s) is 
a function of a particular value of the deterministic parameter vector 0. 

"Note that Jc(<?i) does not depend on any particular value of 61 for a 
random parameter vector 0i, whereas it is a function of 0i for a deterministic 
parameter vector d\. 

'''The light hand side of )13t is known as the Schur complement of the 
matrix C |66|. 

"When 81 £ R^, the EFIM has only one element since ,Js{9\) is a scalai'. 



SHEN AND WIN: FUNDAMENTAL LIMITS OF WIDEBAND LOCALIZATION— PART I: A GENERAL FRAMEWORK 



5 



III. Evaluation of EFIM 



In this section, we apply the notion of equivalent Fisher 
information (EFI) to derive the SPEB for both the case with 
and without a priori knowledge. We also introduce the notion 
of ranging information (RI), which turns out to be the basic 
component of the SPEB. 

A. EFIM without A Priori Knowledge 

First consider a case in which a priori knowledge is unavail- 
able. We apply the notion of EFI to reduce the dimension of 
the original FIM in and the EFIM for the agent's position 
is presented in the following proposition. 

Proposition 1: When a priori knowledge is unavailable, an 
EFIM for the agent's position is 

Je(p,{/^fc : fceA^L}) = 4 TlAlTJ, (14) 



where Tl and Al are given by (ITTT i and ( |42] |. respectively. 

Proof: Let A = TnlAnlTSl + TlAlTT B = TnlAnl, 
and C = Anl in ©■ Applying the notion of EFI in ( fT3] ) leads 
to the result. □ 

Remark 1: When a priori knowledge is unavailable, NLOS 
signals do not contribute to the EFIM for the agent's position. 
Hence we can eliminate these NLOS signals when analyzing 
localization accuracy. This observation agrees with the results 
of ||29l , but the amplitudes of the MFCs are assumed to be 
known in their model. 

Note that the dimension of the EFIM in ( fT4l i is much 
larger than 2x2. We will apply the notion of EFI again to 
further reduce the dimension of the EFIM in the following 
theorem. Before the theorem, we introduce the notion of the 
first contiguous-cluster and ranging information (RI). 

Definition 3 (First Contiguous-Cluster): The first 
contiguous-cluster is defined to be the set of paths 
{1, 2, . . . , j}, such that |tj -Tj+i| < Tj for z = 1, 2, . . . , j - 1, 
and \Tj — Tj+i| > Ts, where is the duration of s{t). 

Definition 4 (Ranging Information): The ranging informa- 
tion (RI) is a 2 X 2 matrix of the form A Jr(0), where A is a 
nonnegative number called the ranging information intensity 
(RII), and Jr(0) a 2 x 2 matrix called the ranging direction 
matrix (RDM) with angle 0, given by 



J.-(0) = 



cos (f) 

cos 6 sin ( 



cos <p sm ( 

,,:„2 / 



The first contiguous-cluster is the first group of non-disjoint 
paths (see Fig. |2]|0 The RDM is one-dimensional along 
the direction with unit intensity, i.e., Jr('/') has one (and 
only one) non-zero eigenvalue equal to 1 with corresponding 
eigenvector q = [ cos (f) sin </> ] . 

Theorem 1: When a priori knowledge is unavailable, the 
EFIM for the agent's position is a 2 x 2 matrix 



Je(p) 



fcGA/l 



Afe Jr {(j)k) , 



(15) 



'^The first contiguous-cluster, defined for general wideband received sig- 
nals, may contain many MFCs. Two paths that arrive at time and Tj are 
called non-disjointed if \Ti — Tj\ < Ts. 




not overlap 



first contiguous-cluster 



Fig. 2. An illustration of the first contiguous-cluster (containing I paths) in 
a LOS signal. 



where A^ is the RII from anchor k, given by 



Afe — 



^(l-Xfc)SNR« 



In (O, < Xfe < 1 is given by ( 1611 ). 



1/2 



C^\S{fWdf 



and 



SNR 



(16) 



(17) 



(18) 



Furthermore, only the first contiguous-cluster of LOS signals 
contains information for localization. 

Proof: See Appendix IC-AI □ 
Remark 2: In Theorem [T] /3 is known as the effective 
bandwidth fSOl, f67l, Xk is called path-overlap coefficient 
(POC) that characterizes the effect of multipath propagation 
for locaHzation, and SNR^'^ is the SNR of the khpath in rfe(t). 
We draw the following observations from Theorem [T] 

> The original FIM in (|7]i can be transformed into a simple 
2x2 EFIM in a canonical form, given by ( fTSl ), as a 
weighted sum of the RDM from individual anchors. Each 
anchor (e.g. anchor k) can provide only one-dimensional 
RI along the direction (pk, from the anchor to the agent, 
with intensity A^ 

• The RII Afc depends on the effective bandwidth of s{t), 
the SNR of the first path, and the POC. Since < Xfe < 1, 
path-overlap in the first contiguous-cluster will reduce the 
RII, thus leading to a higher SPEB, unless the signal via 
the first path does not overlap with others (xk = 0). 

• The POC Xk in ( l6TT i is determined only by the waveform 
s(t) and the NLOS biases of the MFCs in the first 
contiguous-cluster The independence of Xfe on the path 
amplitudes seems counter-intuitive. However, this is due 
to the fact that, although large a^'^ causes severe interpath 
interference for estimating the TOA t^^', it increases 

"For notational convenience, we suppress the dependence of 0^ and 
on the agent's position p throughout the paper. 
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the estimation accuracy for r^.'\ which in turn helps to 
mitigate the interpath interference. 
We can speciahze the above theorem into a case in which 
the first path in a LOS signal is completely resolvable, i.e., 
the first contiguous-cluster contains only a single component. 

Corollary 1: When a priori knowledge is unavailable and 
the first contiguous-cluster of the received waveform from 
anchor k contains only the first path, the RII becomes 



(19) 



Proof: See Appendix |C-B| □ 
Remark 3: When the first path is resolvable, = in ( fT6] l 
and hence attains its maximum value. However, when the 
signal via other paths overlap with the first one, these paths 
will degrade the estimation accuracy of the first path's arrival 
time and hence the RII. Corollary [1] is intuitive and important: 
the RII of a LOS signal depends only on the first path if the 
first path is resolvable. In such a case, all other paths can be 
eliminated, and the multipath signal is equivalent to a signal 
with only the first path for localization. 

From Theorem[T] the SPEB can be derived in ( |20l ) shown at 
the bottom of the page. When the first paths are resolvable, by 
Corollary [H we have all Xk = in (l20b and the corresponding 
7'(p) becomes the same as those based on single-path signal 
models in |[9|, ll29l . However, those results are not accurate 
when the first path is not resolvable. 

B. EFIM with A Priori Knowledge 

We now consider the case where there is a priori knowledge 
of the channel parameters, but not of the agent's position. 
In such cases, since p is deterministic but unknown, /(p) 
is eliminated in ([8]). Similar to the analysis in the previous 
section, we can derive the 2x2 EFIM for the corresponding 
FIM in ([Toll. 

Theorem 2: When a priori knowledge of the channel pa- 
rameters is available and the sets of channel parameters 
corresponding to different anchors are mutually independent, 
the EFIM for the agent's position is a 2 x 2 matrix 

Je(p)= AfcJi-(0fc)+ ^k:iA4>k) , (21) 

feeA^L fceA/'NL 

where is given by ( |65a| ) for LOS signals and ( |65bt for 
NLOS signals. 

Proof: See Appendix FC-CI □ 
Remark 4: Theorem |2] generalizes the result of Theorem [T] 
from deterministic to hybrid parameter estimation^ In this 
case, the EFIM can still be expressed in a canonical form as 
a weighed sum of the RDMs from individual anchors. Note 
that due to the existence of a priori channel knowledge, the 

'^This is the case where the agent's position p is deterministic and the 
channel parameters are random. 



RII of NLOS signals can be positive, and hence these signals 
contribute to the EFIM as opposed to the case in Theorem [T] 

Corollary 2: A priori channel knowledge increases the RII. 
In the absence of such knowledge, the expressions of RII in 
( 165 al l and ( I65bb reduce to (fT6l l and zero, respectively. 

Proof: See Appendix fC-DI □ 

Corollary 3: LOS signals can be treated as NLOS signals 
with infinite a priori Fisher information of b^i}\ i.e., 6^^^ 



is known. Mathematically, ( 165 ab is equivalent to i65H with 
Fe{e;bi'\bi}^)^oo. 

Proof: See Appendix IC-EI □ 

Remark 5: Corollary |2] shows that Theorem |2] degenerates 
to Theorem [T] when a priori channel knowledge is unavailable. 
Moreover, Corollary |3] unifies the LOS and NLOS signals 
under the Bayesian estimation framework: the NLOS biases 
fe^^"* (fc G A/l) can be regarded as random parameters with 
infinite a priori Fisher information instead of being eliminated 
from 9 as in Section III-AI Hence, all of the signals can be 
modeled as NLOS, and infinite a priori Fisher information of 
b\^^ will be assigned for LOS signals. 

We next consider the case where a priori knowledge of the 
agent's position is available in addition to channel parameters. 
Note that the topology of the anchors and the agent changes 
with the position of the agent. The 2x2 EFIM is given in 
(|67] |. which is more intricate than the previous two cases. To 
gain some insights, we consider a special case wherJ^ 



^p{5(p)} = .9(P) , 



(22) 



in which p = Ep {p} is the agent's expected position, for 
some function g{-) involved in the derivation of the EFIM 
(see Appendix IC-Fb . 

Proposition 2: When the a priori position distribution of 
the agent satisfies ( |22] |. and the sets of channel parameters 
corresponding to different anchors are mutually independent, 
the EFIM for the agent's position is a 2 x 2 matrix 



Je(p) 



E 



; Jr {4>k) + Sp 



(23) 



where A^ is given by ( l68b . and (f)k is the angle from anchor k 
to p. 

Proof: See Appendix IC-FI □ 
Remark 6: The a priori knowledge of the agent's position 
is exploited, in addition to that of the channel parameters, for 
localization in Proposition |2] The expressions for the EFIM 
can be involved in general. Fortunately, if (|22] | is satisfied, 
the EFIM can be simply written as the sum of two parts as 
shown in ( 1231 ): a weighted sum of the RDMs from individual 
anchors as in the previous two cases, and the EFIM from 

"This occurs when the agent's a priori position distribution is concentrated 
in a small area relative to the distance between the agent and the anchors, 
so that g(p) is flat in that area. For example, this condition is satisfied in 
far-field scenarios. 



2E,eA^,(l 



^'^'^^ EkeJ^, E,„e^,(l - Xk){l - Xm) SNrWsNRW sin^ 



(20) 
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the a priori knowledge of the agent's position. This resuh 
unifies the contribution from anchors and that from the a priori 
knowledge of the agent's position into the EFIM. The concept 
of localization with a priori knowledge of the agent's position 
is useful for a wide range of applications such as successive 
localization or tracking. 

IV. Wideband Localization with Antenna Arrays 

In this section, we consider localization systems using 
wideband antenna arrays, which can perform both TOA and 
AOA measurements. Since the orientation of the array may be 
unknown, we propose a model to jointly estimate the agent's 
position and orientation, and derive the SPEB and the squared 
orientation error bound (SOEB). 

A. System Model and Squared Orientation Error Bound 

Consider a network where each agent is equipped with an 
Ai"a-antenna array0 which can extract both the TOA and AOA 
information with respect to neighboring anchors. Let A/l — 
{1, 2, . . . , A/'a} denote the set of antennas, and let pn"^^ = 
[a;™™'' y^'^'^Y denote the position of the agent's nth antenna, 
which needs to be estimated. Let (pnk denote the angle from 
anchor k to the n knth antenna, i.e.. 



tan 



-1 Vn 



AiTay 



Vk 



X 



AiTay 



Since relative positions of the antennas in the array are usually 
known, if we denote p = [a; 
as the orientation of the array 
antenna in the array can be represented as (Fig. O 



^ as a reference point and tp 
then the position of the nth 



A2/„(p,V3) 



e A/; 



where Ax„(p, ip) and Aj/„(p, ip) denote the relative distance 
in X and y direction from the reference point to the nth 
antenna, respectively. 

Since the array orientation may be unknown, we classify the 
localization problem into orientation-aware and orientation- 
unaware cases, where ip can be thought of as a random pa- 
rameter with infinite (orientation-aware) and zero (orientation- 
unaware) a priori Fisher information |46|. 

The received waveform at the agent's nth antenna from 
anchor k can be written as 



r„kit) 



where a^^l, and r^'^ are the amplitude and delay, respectively, 
of the kh path, Lnk is the number of MFCs, and Znk{t) repre- 
sents the observation noise modeled as additive white Gaussian 
processes with two-side power spectral density The 



1=1 



.(0 



+ Znk{t), t6[0,Tob), 



Each anchor ha.s only one antenna here. We will di.scu.ss the case of 
multiple antennas on anchors at the end of this section. 

-'Note from geometry that the orientation ifi is independent of the specific 
reference point. 



AiTay 

Pa 



p Array 



Array 

Pi 



\ Array 

Pn. 



Fig. 3. An antenna array is described by the reference point p, the orientation 
ifi, and the relative positions of the antennas. 



relationship between the position of the nth antenna and the 
delay of the kh path is 
1 



.(0 



AiTay 



Pfel 



u(0 



(24) 



The parameters to be considered include the position of 
the reference point, the array orientation, and the nuisance 
multipath parameter as 

- T - T 1 T 



e 



"N., 



(25) 



where /«„ consists of the multipath parameters associated with 
the received waveforms from all anchors at the nth antenna, 

T T 1 T 



and each Knk consists of the multipath parameters associated 
with r„fe(t). 



^nk 



..(1) 



nk 



a, 



nk 



Similar to Section III-BI the overall received waveforms at 
the antenna array can be represented, using the KL expansion, 
by r = [ r| r| • • • ] , where 

~ [ r^,l ^Ji,2 ■ ■ • ^Ji.Nb ] 7 

in which r„fe is obtained by the KL expansion of rnk{t)- 
Definition 5 ( Squared Orientation Error Bound): The 
squared orientation error bound (SOEB) is defined to be 

B. EFIM without A Priori Knowledge 

We first consider scenarios in which a priori knowledge is 
unavailable. Following similar steps in Section IIII-BI we have 
the following theorem. 

Theorem 3: When a priori knowledge is unavailable, the 
EFIMs for the position and the orientation, using an A^^- 
antenna array, are given respectively by 



jA„-ay(p) ^ ^ j^(pA.ay) _ ^ 1 

nGA/'a 2^neM, 2^ki 



qq 



(26) 
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and 



(27) 



where A„fe is given by d?!]), q„fc = [cost 
and 

q = E E '^"'^ ' 

and 



(28) 



hnk = X" ^2^«(P, COS (fink + -J- Aj/„(p, (p) Sin0„fc 



Proof: See Appendix ID-AI 
Corollary 4: The EFIM for the position is given by 



(29) 

□ 

(30) 



for orientation-aware locaHzation. 

Proof: (OutHne) In orientation-aware locaHzation, the 
angle ip is known and hence excluded from the parameter 
vector 6 in dZST l. Consequently, the proof of this corollary is 
analogous to that of Theorem [3] except that the components 
corresponding to ip are eliminated from the FIM in i69[ and 
(ITOI i. One can obtain (|30] | after some algebra. □ 
Remark 7: The EFIM Je(p^™'') in ^ and ^ coiTe- 
sponds to the localization information from the nth antenna. 
We draw the following observation from the above theorem. 
• The EFIM J^"^^(p) in (l26T l consists of two parts: 1) 
the sum of localization information obtained by indi- 
vidual antennas, and 2) the information reduction due 
to the uncertainty in the orientation estimate, which 
is subtracted from the first part@ Since qq^ in the 
second part is a positive semi-definite 2x2 matrix and 
E„GA/-, EfceA/-, ^nfc^L is always positive, we have the 
following inequality 



jf"^^(p)^ E JeCpr")- 



(31) 



The inequality implies that the EFIM for the position, 
using antenna arrays, is bounded above by the sum of 
all EFIMs corresponding to individual antennas, since 
the uncertainty in the orientation estimate degrades the 
localization accuracy, except for q = or orientation- 
aware localization (i.e., dSOll). 
• The EFIM J^"^^(p) and J(f™''(<^) depend only on the 
individual RI between each pair of anchors and anten- 
nas (through Anfc's and (/)„fc's), and the array geometry 
(through /infe's). Hence it is not necessary to jointly 
consider the received waveforms at the antennas, 
implying that AOA obtained by antenna arrays does not 

^^For notational convenience, we suppress the dependence of A„j., 
and q on the reference position p. 



increase position accuracy. Though counterintuitive at 
first, this finding should not be too surprising since AOA 
is obtained indirectly by the antenna array through TOA 
measurements, whereas the TOA information has already 
been fully utilized for localization by individual antennas. 
• The gain of using antenna arrays for localization mainly 
comes from the multiple copies of the waveform received 
at the A'a antennas (see (E^ll)!^ and its performance is 
similar to that of a single antenna with A^a measurements. 
The advantage of using antenna arrays lies in their ability 
of simultaneous measurements at the agent. 

The equality in dSTT i is always achieved, independent of ref- 
erence point, in orientation-aware localization. However, only 
a unique reference point achieves this equality in orientation- 
unaware localization. We define this unique point as the 
orientation center 

Definition 6 ( Orientation Center): The orientation center is 
a reference point p* such that 

jAn-ay(p*)^ E ^Pn'"') ' 

Proposition 3: Orientation center p* exists and is unique in 
orientation-unaware localization, and hence for any p 7^ p*, 



Proof: See Appendix fP-BI □ 
Remark 8: The orientation center p* generally depends on 
the topology of the anchors and the agent, the properties of 
the received waveforms, the array geometry, and the array 
orientation. Since q = at the orientation center, the EFIMs 
for the array center and the orientation do not depend on 
each other, and hence the SPEB and SOEB can be calculated 
separately. The proposition also implies that the SPEB of 
reference points other than p* will be strictly larger than that 
of p*. The SPEB for any reference point is given in the next 
theorem. 

Corollary 5: The SOEB V{(p) is independent of the refer- 
ence point p, and the SPEB is 



nP)=^(P*) + i|p-P*f •T'M. 



(32) 



Proof: See Appendix fP-CI □ 
Remark 9: The SOEB does not depend on the specific 
reference point, which was not apparent in (l27b . However, 
this is intuitive since different reference points only introduce 
different translations, but not rotations. On the other hand, 
different reference point p results in different /i„fe's and hence 
different q, which in turn gives different EFIM for position 
(see (|26]l). We can interpret the relationship in ( [32] l as follows: 
the SPEB of reference point p is equal to that of the orientation 
center p* plus the orientation-induced position error, which is 
proportional to both the squared distance from p to p* and 
the SOEB. 

^^In near-field scenarios where the antenna separation is on the order of 
the distances between the array and the anchors, additional gain that arises 
from the spatial diversity of the multiple antennas may be possible. 
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C. EFIM with A Priori Knowledge 

We now consider a scenario in which the channel parameter 
vector Knk independent for different n's and fc's. The inde- 
pendence assumption serve as a reasonable approximation of 
many realistic scenarios, especially near-field cases. When the 
different sets of channel parameters are correlated, our results 
provide an upper bound for the EFIM. 

Proposition 4: When a priori knowledge of channel pa- 
rameters is available and the set of channel parameters cor- 
responding to different anchors and antennas are mutually 
independent, the RII Xnk becomes ( |72| |. 

Proof: See Appendix ID-AI □ 

We then consider the case where a priori knowledge of 
the agent's position and orientation is available in addition 
to channel knowledge. Note that the topology of the agent's 
antennas and anchors changes with the agent's positions and 
orientations. The expression of the EFIM can be derived 
analogous to ( IFtT i. which is involved in general. Again to 
gain insights about the contribution of a priori position and 
orientation knowledge, we consider scenarios under condition 

^P,v{9{P,'P)} = 9{P:'f) , (33) 

where <p = E,^ {ip}, for some functions g{-) involved in the 
derivation of the EFIM. 

Corollary 6: When a priori position and orientation dis- 
tribution of the agent satisfies (|33] |. and the sets of channel 
parameters corresponding to different anchors and antennas 
are mutually independent, the EFIMs for the position and the 
orientation, using an A'a-antenna array, are given respectively 
by 

jf-^(p)=Hp+ J2 E '^nkJMnk) 

neAf-i keAfb 



and 

E "Y >'nkJr{4>nk) +ap\ q, 

VneM, ke^fb ) 

where A„fc, 0„fc, hnk and q are corresponding functions in 
Theorem [3] of p and (p, respectively, and = Y'e(Lp\ tp^ ip). 

Proof: (Outline) The proof of this corollary is analogous 
to that of Theorem [3] Note that when condition ( |33] | is 
satisfied, the a priori knowledge of position and orientation 
for localization can be characterized in the EFIM by using the 
approximation as in the proof of Proposition |2] □ 

D. Discussions 

1 ) Far-field scenarios: The antennas in the array are closely 
located in far-field scenarios, such that the received waveforms 
from each anchor experience statistically the same propagation 
channels. Hence we have (j)nk = (f'k and A„fc — Xk for all 
n, leading to Je(pn"^'^) = Je(p)- We define an important 
reference point as follows. 



Definition 7 (Array Center): The array center is defined as 
the position po, satisfying 

^ A2;„(po, <p) = and ^ A?/„(po, (p) = . 

Proposition 5: The array center becomes the orientation 
center in far-field scenarios. 

Proof: See Appendix ID^ □ 

Remark 10: Since the orientation center has the minimum 
SPEB, Proposition |5] implies that the array center always 
achieves the minimum SPEB in far-field scenarios. Hence 
the array center is a well-suited choice for the reference 
point, since its position can be determined from the array 
geometry alone, without requiring the received waveforms and 
the knowledge of the anchor's topology. 

In far-field scenarios, we choose the array center po as the 
reference point. The results of Theorem [3] become 

Jr^''(P0)=iVa E ^kU^k), 

keN'b 

and 

J^"''i^)^ E E ^^~^nk, 

where h„k is a function of po. Similarly, when the a priori 
position and orientation knowledge is available and condition 
(|33] | is satisfied, the results of Proposition |6] become 

J^^^(po)=7Va E AfcJr(0fc)+Hp, 

ke^^b 

and 

Jt'^'i^)- E E^fc^«fc+=^' 

where hnk is a function of po = Epi-,{po}. 

Note that the localization performance of an A^a-antenna 
array is equivalent to that of a single antenna with 
measurements, regardless of the array geometry, in far-field 
scenarios. 

2) Multiple antennas at anchors: When anchors are 
equipped with multiple antennas, each antenna can be viewed 
as an individual anchor. In this case, the agent's SPEB goes 
down with the number of the antennas at each anchor. Note 
that all the antennas of a given anchor provide RI approxi- 
mately in the same direction with the same intensity, as they 
are closely located. 

3) Other related issues: Other issues related to localization 
using wideband antenna arrays include the AOA estimation, 
the effect of multipath geometry, and the effect of array 
geometries. A more comprehensive performance analysis can 
be found in LI IJ . 

V. Effect of Clock Asynchronism 

In this section, we consider scenarios in which the clocks of 
all anchors are perfectly synchronized but the agent operates 
asynchronously with the anchors ll68l . In such a scenario, the 
one-way time-of-flight measurement contains a time offset 
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between the agent's clock and the anchors' clock0 Here, 
we investigate the effect of the time offset on localization 
accuracy. 

A. Localization with a Single Antenna 

Consider the scenario described in Section [III where each 
agent is equipped with a single antenna. When the agent 
operates asynchronously with the anchors, the relationship of 
(|2|i becomes 

1 



Pfcl 



+ B 



where S is a random parameter that characterizes the time 
offset in terms of distance, and the corresponding parameter 
vector 6 becomes 



6/ = [ pT B 



Similar to Theorem |2] where p is deterministic but unknown 
and the remaining parameters are random, we have the fol- 
lowing result. 

Theorem 4: When a priori knowledge of the channel pa- 
rameters and the time offset is available, and the sets of 
channel parameters corresponding to different anchors are 
mutually independent, the EFIMs for the position and the time 
offset are given respectively by 



(34) 



and 



UB) = Afc + Sb - eg 




qB , (35) 



where Afc is given by ( I65bl i. Qb J2keAfb 

Sb = F£)(i? ; B, B) . 

Proof: See Appendix lE-AI □ 
Remark 11: Since Qb Qb is a positive semi-definite matrix 
and X]feeA4 positive in ( l34l l. compare to Theorem |2] we 
always have the inequality 

Je(P)^Je(p), (36) 

where the equality in (l36T l is achieved for time-offset-known 
localization (i.e., Sb — oo), or time-offset-independent local- 
ization (i.e., Qb — 0). The former corresponds to the case 
where accurate knowledge of the time offset is available, 
while the latter depends on the RII from each anchor, and 
the topology of the anchors and agent. The inequality of (|36] | 

^''We consider scenarios in wliicli localization time is short relative to clock 
drifts, such that the time offset is the same for all measurements from the 
anchors. 



results from the uncertainty in the additional parameter B, 
which degrades the localization accuracy. Hence the SPEB in 
the presence of uncertain time offset is always larger than or 
equal to that without a offset or with a known offset. 

We next consider the case where a priori knowledge of 
the agent's position is available. When the a priori position 
distribution of the agent satisfies (|22] |. we have the following 
corollary. 

Corollary 7: When the a priori position distribution of 
the agent satisfies (l22l l. and the sets of channel parameters 
corresponding to different anchors are mutually independent, 
the EFIMs for the position and the time offset are given 
respectively by 



E 



— qe Qb 



keMh 



and 



UB) E ^fc + ^B-q^ (Hp+ E 



Afc +S 



Afe Jr(0fe) qB , 



where 4>k is the angle from anchor k to p, Afc is given by ([68]), 
and cjb is a function of p. 

Proof: (Outline) Conditions in ( l22b hold in far-field sce- 
narios, and we can approximate the expectation over random 
parameter vector p using the average position p. By following 
the steps of Theorem 2] and Proposition |2] we can derive the 
theorem after some algebra. □ 

B. Localization with Antenna Arrays 

Consider the scenario describing in Section |IV] where each 
agent is equipped with an array of N.^ antennas. Incorporating 
the time offset B, (|24] | becomes 

1 



^(0 

'nk 



^Anay 



Pfcl 



^nk 



B 



and the corresponding parameter vector 6 becomes 

6 — [ ip B hi[ ■ ■ ■ kjf ] ^ . 

Similar to Theorem [3] where p and ip are deterministic but 
unknown and the remaining parameters are random, we have 
the following theorem. 

Theorem 5: When a priori knowledge of the channel pa- 
rameters is available, and the sets of channel parameters 
corresponding to different anchors and antennas are mutually 
independent, the EFIM for the position, the orientation, and 
the time offset, using an A^a-antenna array, is given by dJTl ) 
shown at the bottom of the page, where S = cx) and S = 
correspond to orientation-aware and orientation-unaware lo- 
calization, respectively, and \nk, ^.nk, and hnk are given by 
( |72l ), ( |28T l, and respectively. 



E n eK, EfcGA^ ^nk Kik qnfc E n 
•^^ ^ = EnSA/;, EfceAAb ^nk \ik q^fe EneX EfceA^ ^"nk ^nk + ^ EnGA/;, EfceA4 ^^k ^nk 



(37) 
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Proof: See Appendix lE-BI □ 
Remark 12: Theorem|5]gives the overall 4x4 EFIM for the 
position, the orientation, and the time offset, where individual 
EFIMs can be derived by applying the notion of EFI again. 

We finally consider the case where a priori knowledge of 
the agent's position and orientation is available. The EFIM in 
far-field scenarios is given in the following corollary. 

Corollary 8: When a priori knowledge of the agent's po- 
sition, orientation, time-offset, and the channel parameters is 
available, and the sets of channel parameters corresponding to 
different anchors and antennas are mutually independent, in 
far-field scenarios, the EFIMs for the position, the orientation, 
and the time offset, using an A^a-antenna array, are given 
respectively by 

J«(po)-iV. E AfeJr(^A0+3p 

fceA/b 

1 



^A.ay-B(^)^ E E 



and 



-Qb f^a ( 



Qb 



where po is the expected position of the agent's array center, 
is the angle from anchor fc to po, and A^, Qb, and are 
functions of po. 

Proof: See Appendix lE-CI □ 
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Fig. 4. Network topology: Four anchors are equally spaced on a circle with 
an agent at the center. All signals from the anchors to the agent are LOS. 



a priori knowledge f50l. High SNR can be attained using se- 
quences with good correlation properties |69|-|71|, or simply 
repeated transmissions. Therefore, the SPEB is achievable. 

C. Generalization to 3D Localization 

All results obtained thus far can be easily extended to three- 
dimensional case, i.e., p = [x y and the RDM becomes 

where ipk and (pk are the angles in the coordinates, and 



qfc 



cos (fik cos cpk sm ifk cos (pk sm cpk 



Similarly, we can obtain a corresponding 3x3 EFIM in the 
form of (EB. 



VI. Discussions 

In this section, we will provide discussions on some related 
issues in the paper It includes 1) the relations of our results to 
the bounds based on signal metrics, 2) the achievability of the 
SPEB, and 3) generalization of the results to 3D localization. 

A. Relation to Bounds Based on Signal Metrics 

Analysis of localization performance in the literature mainly 
employs specific signal metrics, such as TOA, AOA, RSS, and 
TDOA, rather than utilizing the entire received waveforms. 
Our analysis is based on the received waveforms and exploits 
all the localization information inherent in these signal metrics, 
implicitly or explicitly. In particular, TOA and joint TO A/ AOA 
metrics were incorporated in our analysis in SectionlllllandllVI 
respectively. Similarly, TDOA and joint TDOA/ AOA metrics 
were included in the analysis of Section[Vl and the RSS metric 
has been implicitly exploited from a priori channel knowledge 
in Section lirC2l 

B. Achievability of the SPEB 

Maximum a posterior (MAP) and maximum likelihood 
(ML) estimates respectively achieve the CRB asymptotically 
in the high SNR regimes for both the case with and without 



VII. Numerical Results 

In this section, we illustrate applications of our analytical 
results using numerical examples. We deliberately restrict our 
attention to a simple network to gain insights, although our 
analytical results are valid for arbitrary topology with any 
number of anchors and any number of MPCs in the received 
waveforms. 

A. Effect of Path-Overlap 

We first investigate the effect of path-overlap on the SPEB 
when a priori knowledge is unavailable. In particular, we 
compare the SPEB obtained by the full-parameter model 
proposed in this paper and that obtained by the partial- 
parameter model proposed in |281. In the partial-parameter 
model, the amplitudes of MPCs are assumed to be known and 
hence excluded from the parameter vector. 

Consider a simple network with four anchors (iVb — 4) 
equally spaced on a circle and an agent at the center receiving 
all LOS signals (see Fig. |4|i. Each waveform consists of two 
paths: one LOS path (SNr!.^' = dB) and one NLOS path 

(2) 

(SNR^ = —3 dB), and the separations of the two paths 
T^^^ — T^/'' are identical for all k. In addition, the transmitted 
waveform is a second derivative of Gaussian pulse with width 
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Full-parameter 

Partial-parameter 

Non-overlap 




0.01 



1 2 3 4 5 

Path separation (ns) 

Fig. 5. SPEB as a function of path separation for the full-parameter, partial- 
parameter, and non-overlap models, without a priori knowledge. 



approximately equal to 4 ns. Figure |5] shows the SPEB as a 
function of path separation r^^^ — rj,^^ according to Theorem 

m 

We can draw the following observations. First of all, path- 
overlap increases the SPEB in both models, since it reduces the 
ability to estimate the first path and hence decreases the RII. 
Note that the shape of the curves depends on the autocorrela- 
tion function of the waveform s{t) |47|. Secondly, when the 
path separation exceeds the pulse width (approximately 4 ns), 
the two models give the same SPEB, which equals the non- 
overlapping case. In such cases, the first contiguous-cluster 
contains only the first path, and hence the RII is determined 
by the first paths. This agrees with the analysis in Section Hill 
Thirdly, excluding the amplitudes from the parameter vector 
incorrectly provides more RI when the two paths overlap, and 
hence the partial-parameter model results in a loose bound. 
This demonstrates the importance of using the full-parameter 
model. 



1 2 3 4 5 

Path separation (ns) 

(a) SPEB with a priori knowledge of q^,^' and while H (2) = 0. 




1 2 3 4 5 

Path separation (ns) 



(b) SPEB with a priori knowledge of b, , while H (d = S (21 = 0. 

Fig. 6. SPEB with a priori iaiowledge of the amplitudes and the NLOS 
biases as a function of path separation, respectively. 



B. Improvement from A Priori Channel Knowledge 

We then quantify the contribution of the a priori knowledge 
of channel parameters to the SPEB. The network topology 
and channel parameters are the same as those in Section 
IVII-AI except a priori knowledge of ol]^\ Ojf' and fo^'^'' is now 
available. For simplicity, we consider these parameters to be 
independent a priori and denote the a priori Fisher information 
of parameter 9i by Sg^ — Fe{9 ;9i,9i). In Fig. |6(a)[ the 
SPEBs are plotted as functions of the path separation for 
different a priori knowledge of a[.^'' and a[,^^ (no a priori 
knowledge of fo^^''); while in Fig. |6(b)| the SPEBs are plotted 
for different a priori knowledge of b^, ' (no a priori knowledge 
of a^^^ and a^^^). 

We have the following observations. First of all, the SPEB 
decreases with the a priori knowledge of the amplitudes and 
the NLOS biases. This should be expected since a priori 
channel knowledge increases the RII and thus localization 



accuracy, as indicated in Corollary |2] Moreover, the NLOS 
components are shown to be beneficial for localization in the 
presence of a priori biases knowledge, as proven in Section 
IIII-BI Secondly, as the a priori knowledge of the amplitudes 
approaches infinity, the SPEB in Fig. |6(a)| obtained using the 
full-parameter model converges to that in Fig. |5] obtained 
using the partial-parameter model. This is because the partial- 
parameter model excludes the amplitudes from the parameter 
vector, which is equivalent to assuming known amplitudes and 
hence infinite a priori Fisher information for the amplitudes 
(S (1) = S (2) = 00). Thirdly, it is surprising to observe that, 
when the a priori knowledge of the NLOS biases is available, 
path-overlap can result in a lower SPEB compared to non- 
overlapping scenarios. This occurs at certain regions of path 
separations, depending on the autocorrelation function of s{t). 
Intuitively, path-overlap can lead to a higher SNR compared to 
non-overlapping cases, when a priori knowledge of the NLOS 
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0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
Path inter-arrival time 1/v (ns) 



Fig. 7. POC as a function of the path inter-arrival time for different number 
of MPCs. 



biases is available. 

C. Path-Overlap Coefficient 

We now investigate the dependence of POC x on path 
arrival rate. We first generate channels with L MPCs according 
to a simple Poisson model with a fixed arrival rate i/, and then 
calculate x according to i6l\ . Figure |7] shows the average path- 
overlap coefficient as a function of path inter-arrival rate (l/i^) 
for different number of MPCs, where the averaging is obtained 
by Monte-Carlo simulations. 

We have the following observations. First of all, the POC x 
is monotonically decreasing from 1 to with l/iy. This agrees 
with our intuition that denser multipath propagation causes 
more interference between the first path and other MPCs, and 
hence the received waveform provides less RII. Secondly, for 
a fixed i/, the POC increases with L. This should be expected 
as additional MPCs may interfere with earlier paths, which 
degrades the estimation accuracy of the first path and thus 
reduces the RII. Thirdly, observe that beyond L = 5 paths, x 
does not increase significantly. This indicates that the effect of 
additional MPCs beyond the fifth path do not contribute to the 
RII, regardless of the power dispersion profile of the received 
waveforms. 

D. Outage in Ranging Ability 

We have observed that the channel quality for ranging is 
characterized by the POC. If the multipath propagation has a 
larger POC (close to 1), we may consider the channel in outage 
for ranging. We define the ranging ability outage (RAO) as 

Pout(xth) = P{x > Xth}, 

where Xth is the threshold for the POC. The RAO tells us 
that with probability Pout(Xth), the propagation channel is 
unsatisfactory for ranging. 

The RAO as a function of Xth for different Poisson arrival 
rate is plotted in Fig. [8] for a channel with L = 50. The RAO 




0.2 0.4 0.6 0.8 1 

Threshold xih 



Fig. 8. RAO as a function of the threshold Xth for different path inter- 
arrival time 1/v wfith L = 50. The five curves correspond to inter-arrival 
time 1/v = 3.5, 2.5, 2, 1.6, 1.4 ns, respectively. 



decreases from 1 to 0, as the threshold Xth increases or the 
path arrival rate v decreases. This should be expected because 
the probability of path-overlap decreases with the path arrival 
rate, and consequently decreases the RAO. The RAO can be 
used as a measure to quantify the channel quality for ranging 
and to guide the design of the optimal transmitted waveform 
for ranging. 

E. SPEB and SOEB for Wideband Antenna Array Systems 

We consider the SPEB and SOEB for different reference 
points of a uniform linear array (ULA). The numerical results 
are based on a network with six equally spaced anchor nodes 
(A^b = 6) located on a circle with an agent in the center The 
agent is equipped with a 4-antenna array {N-^ — 4) whose 
spacing is 0.5 m. In far-field scenarios, A„fc = = 10 and 
</>„fe = (j)f,. Fig. |9(a)]and Fig. |9(b)] show the SPEB and SOEB, 
respectively, as a function of different reference point along 
the ULA for different a priori knowledge of the orientation 
and reference point. 

We have the following observations. First of all, a priori 
knowledge of the orientation improves the localization accu- 
racy as the SPEB decreases with "E.^. The curves for S;^ = 
and = oo correspond to the orientation-unaware and 
orientation-aware cases, respectively. As a counterpart, a priori 
knowledge of the reference point improves the orientation 
accuracy as the SOEB decreases with Hp. This agrees with 
both intuition and Theorem [3] Secondly, the array center has 
the best localization accuracy, and its SPEB does not depend 
on S^, which agrees with Theorem [3] On the other hand, the 
array center exhibits the worst orientation accuracy, and its 
SOEB does not depend on Hp. This should be expected since 
the knowledge for the array center tells nothing about the array 
orientation. Thirdly, the SPEB increases with both the distance 
from the reference point to the array center and the SOEB, as 
predicted by Corollary |5] On the contrary, the SOEB decreases 
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(a) SPEB as a function of the reference point-to-array center distance 



(a) SPEB as a function of Anchor Ai's position 
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Hp=diag{0,0} 
E:p=diag{10,10} 
E:p=diag{l0M02} 
Sp = diag{oC',oo} 
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(b) SOEB as a function of the reference point-to-an'ay center distance 
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No time offset 
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Anchor Ai 's position on the circle 0i (rad) 

(b) STEB as a function of Anchor Ai's position 



Fig. 9. SPEB and SOEB with different a priori Imowledge of agent's position Fig. 10. SPEB and STEB with different a priori Imowledge of the time 



and orientation, respectively 



offset, and = 0, 10, 10 , oo respectively. 



as a function of the distance from the reference point to the 
array center if a priori knowledge of the reference point is 
available. This observation can be verified by Theorem[3] Last 
but not least, the SPEB is independent of specific reference 
point if = oo, as referred to orientation-aware localization, 
and the SOEB is independent of the specific reference point 
if Sp = 0, as shown in Corollary |5] 

F. SPEB with Time Ojfset and Squared Timing Error Bound 

We finally investigate the effect of time offset on the SPEB 
and squared timing error bound (STEB) for the network 
illustrated in Fig. H) The RII from each anchor — 10, 
k G {1,2,3,4}. Initially, four anchors are placed at = 0, 
02 = 7r/2, 03 = TT, and 04 = 37r/2, respectively. We then 
vary the position of Anchor Ai counter-clockwise along the 
circle. Figure [TO(a)| and [TO(b)| show the SPEB and STEB, 
respectively, as functions of 0i for different a priori knowledge 
of the time offset. 



We have the following observations. First of all, both the 
SPEB and STEB decrease with the a priori knowledge of the 
time offset. The SPEB for the case 2b = oo in Fig. |10(a)| 
i.e., known time offset, is equal to that of a system without 
a time offset. On the other hand, when Sb = oo, the STEB 
in Fig. |10(b)| is equal to zero regardless of 0i since the offset 
is completely known. Secondly, all the curves in Fig. |10(a)| 
have the same value at 0i = 0. The time offset has no effect 
on the SPEB at this point, since Qb = 0, referred to as time- 
ojfset- independent localization. In this case, both the SPEB 
and STEB achieve their minimum, implying that location and 
timing information of a network are closely related. Third, 
as 01 increases from to tt, all the curves in Fig. 1 10(a) first 
increase and then decrease, whereas all the curves in Fig. 10(b)| 
increase monotonically. We give the following interpretations: 
the estimation error of time offset in Fig. |10(b)| becomes larger 
when all the anchors tend to gather on one side of the agent (0i 
increases from to tt). In Fig. |10(a)[ the SPEB first increases 
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and 



since both the locahzation information X]fceA4 "^r(</'fe) in as follows 
( [34] l and the information for the time offset becomes smaller. 
Then the SPEB decreases since the localization information 
increases (when > 7r/2) faster compared to the decrease of 
the information for time offset. Note in Fig. |10(a)| that although 
ipi ~ and ipi ~ TT result in the same SPEB in the absence of 
time offset, (pi ^ gives a better performance in the presence 
of time offset. 

VIII. Conclusion 

In this paper, we developed a framework to study wideband 
wireless location-aware networks and determined their local- 
ization accuracy. In particular, we characterized the localiza- 
tion accuracy in terms of a performance measure called the 
squared position error bound (SPEB), and derived the SPEB 
by applying the notion of equivalent Fisher information. This 
methodology provides insights into the essence of the local- 
ization problem by unifying the localization information from 
the a priori knowledge of the agent's position and information 
from individual anchors. We showed that the contributions 
from anchors, incorporating both measurements and a priori 
channel knowledge, can be expressed in a canonical form as a 
weighted sum of the ranging direction matrix. Our results are 
derived from the received waveforms themselves rather than 
the signal metrics extracted from the waveforms. Therefore, 
our framework exploits all the information inherent in the 
received waveforms, and consequently the results in this paper 
serve as fundamental limits of localization accuracy. These 
results can be used as guidelines for localization system 
design, as well as benchmarks for location-aware networks. 

Appendix A 
Fisher Information Matrix Derivation 



-L — 



lnl — 



Gi G2 
Di 

D2 











Da/ 




Al = diag{*i,*2,--- ,*m} 



(41) 



(42) 



Anl = diag{*M+i, *A/+2, • ■ • , *ArJ , 
where = [ hL^^i ], 

Gk=qkll with lfc = [ 1 ■•• 1 ]\ (43) 

2Lfc components 

qA; = [cos(j)k sin(j)k\^, and G M2Lfcx2Lfc given by 



^Fr(r|0;77„r/,). 
Note that elements in '3/^ can be expressed as 
d^lnf{r\e)\ 



(44) 



No 



|27r/^(/)|%xp{-j27r/(r(' 
92 



Si) 



)}df 



No dT^'^drl^^ 



R. 



(4' 



) _ ^0") 



I d"^ In f(r\0) 

To facilitate the analysis, we consider a mapping from E < J \ \ > 



into another parameter vector f] = [ rfl r/l 



where rj/^ = 



4^^ 



n 



T 
T 



with 



= a\. /c. When the agent is localizable^j this mapping 
is a bijection and provides an alternative expression for the 
FIM as 



J j27r/|5(/)|^xp{-j27r/(r«-r(^))}d/ 



2afc d 
No drj:'' 



R. 



(38) 



and 



where J^j is the FIM for t], and T is the Jacobian matrix for 
the transformation from to 77, given respectively by 



92 ln/(r|6») 



J^^Fr(r|0;r7,r7) = 



and 



T ^ ^ = 1 
89 c 



Tl 




Al 


Tnl 
I 





Anl 



(39) 



(40) 



if! 

No 
2c2 

Nr 



l^(/)l'exp{-j27r/(r«-r(^-))}d/ 
where Rs{t) = J s{t)s{t — T)dt. In particular. 



with denoting a matrix of all zeros and I denoting an identity 
matrix. The block matrices Tl, Tnl, Al, and Anl are given 



i^kU - Fr(r|0 ; Ti'\Ti'^) = Sn'fi' SNR^ , 



(46) 



-^Note that an agent is said to be localizable if its position can be 
determined by the signal metrics extracted from the waveforms received from 
neighboring anchors, i.e., triangulation is possible. This is true when M > 3, 
or in some special cases when M = 2. 



where (3 and SNR^*-* are given by (fTTl i and ( fTSl l, respectively. 



Substituting and (gOll into 



, we have the FIM 3g in 
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Appendix B 

Wideband Channel Model and A Priori Channel 
Knowledge 

Wideband channel measurements have shown that MFCs 
follow random arrival and their amplitudes are subject to 
path loss, large and small-scale fading. While our discussion 
is valid for any wideband channels described by (dJ, we 
consider the model of IEEE 802. 15.4a standard for exposition. 
Specifically, this standard uses Poisson arrivals, log-normal 
shadowing, Nakagami small-scale fading with exponential 
power dispersion profile (PDP) [26] . 



A. Path Arrival Time 

The arrival time of MPCs is commonly modeled by a 
Poisson process [|26l . 1641 . Given the path arrival rate v, we 
have 



= V cxp 



for T^'^ > T^'~^^ and I > 2. Using we obtain 



(0 uii-i) 



for b^^> > b 
consistency. 



and / > 1. Note that we let b 



(0) _ 



for 



B. Path Loss and Large-Scale Fading 

The RSS in dB at the distance dk can be written as 

Ffe = Po-10f?logio (^) + w, 



where Pq is the expected RSS at the reference distance 
do, g is the propagation (path gain) exponent, and w is a 
random variable (rv.) that accounts for large-scale fading, or 
shadowing. Shadowing is usually modeled with a log-normal 
distribution, such that w is a Gaussian r.v. with zero-mean and 
vai-iance cr|, i.e., w N (0,cr|)@ The PDF of the RSS of 
rfc(t) can then be written as 



9p{Pk\dk) ocexpi 



Pfc-Po + 10^^1oglo 



(48) 



where dk = ||p ^ PA;||, and Pk is given by 



Pk = 10 log. 



10 



.1=1 



with Es {•} denoting the average over small-scale fading. 

standard deviation is typically 1-2 dB (LOS) and 2-6 dB (NLOS) 
1211 around the path gain. 



C. Power Dispersion Profile and Small-Scale Fading 

As in 1241 . ll26l . we consider an exponential PDP given 



byE!| 



(0 



a,. 



= Qfccxp -- 



7fe 



- QI'^ , (49) 



where •jk is the decay constant, and Qk is a normalization 
coefficient such that 



Qk = 



Ef=oexp(-T^'V7fc) 



(50) 



In addition, af^ is a Nakagami rv. with second moment given 
by (|49T l. Specifically, we have 



ik^dk^Pk 



Mo I D 

mi \ 



exp 



m 
770 



(0 



(51) 

where r(TO/) is the gamma function and mi > 1/2 is the 
Nakagami m-factor, which is a function of ||26l . 

D. A Priori PDF for Multipath Parameters 

The joint PDF of the multipath parameters and the RSS, 
conditioned on the distance from anchor k to the agent, can 
be derived as 

f{ak,hk,Pk \dk) = 9p{Pk\dk) JJ^a* (afe' \hk,dk,Pk 



1=1 



Y[9b. ( 



(52) 



1=1 



By integrating over Pk, we obtain the PDF of the multipath 
parameters of rk{t) as follows 

f{Kk\dk) = f{ak,hk\dk) 

f{cXk,hk,Pk\dk)dPk. (53) 

-OO 



Equation (1531 1 characterizes the a priori knowledge of channel 
parameters, and can be obtained, for IEEE 802.15.4a standard, 
by substituting (|47]l, (08]) and (|5T]) into (|52l) and ^5^. Note that 
since is known, dk is a function of p and hence we have 

Appendix C 
Proofs of the Results in SectionHII] 

A. Proof of Theorem |7] 

Proof: We first prove that Je(p) is given by (flSl l. We 
partition Gfe in ( |43] | and '^k in (l44l i as 

kfc 



Gk = [ qfc Gfc ] and = 



^^Note that the first component of LOS signals can exhibit a stronger 
strength than i49\ in some UWB measurement U2l . In such cases, i49) and 
jSOt need to be modified, accordingly. 
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where [*fc]i,i is obtained by (gS]), e M^Lfc-i^ ^ 
Gfe = qfc [ 1 ••• 1 j"". 

2 L fc — 1 components 

Using these notations, we can write the EFIM given by (fT4l t 
in Proposition [T] after some algebra, in the form of (fT2b , 



+ 5^ {G,kfcqT + q,.kTGT + Gfe*feGl} 
B^[qikT + Gi*i ••• qMkT, + GM*M], 



and 



C^diag{*i,*2,...,*A/} 



Applying the notion of EFI as in ( fT3T l, we obtain the 2x2 
Je(p) as 

87r2/32 



Je(p) = 

where the POC 



^(l-Xfe)SNRWqfcq];, (55) 



Xk 



87r2/32SNR^^^ 



(56) 



This completes the proof of (flST l. 

Next, we show that only the first contiguous-cluster contains 
information for localization. Let us focus on Xfc- Define the 
following notations for convenience: 

Rs{i,i) = Rs{t)\,_H)^(3) , 



and 



d 



If the length of the first contiguous-cluster in the received 
waveform is Lk where I < < Lk, then Rs{i,j) ~ 
Rs{i,j) = Rs{i,j) = for I e {l,2,...,Lfc} and j e 
{Lfe + l,Zfc + 2,...,ifc}, anc@ 



kk^[kl ]^ and *fc ^ 



Kl 



is a block matrix that is irrelevant to the rest of the derivation. 



where k^ e M^Z^-i ^^^j ^ ]^(2L,-i)x(2L,-i)^ q 
becomes 



Xk 



87r2/32SNR^^^ ' 



(57) 



which depends only on the first Lk paths, implying that only 
the first contiguous-cluster of LOS signals contains informa- 
tion for localization. 

Finally, we show that Xk is independent of a^J^K Note that 
ypfc and kfe can be written as 



,c,...,al ',cj Tk 



*fe = ^ diag {c, < 

X diag|c,4^\c,...,af'=\cj, (58) 



and 



(1) 



No 



A- / (2) (L 

diag |c,a;^ ',c,...,al 



tfe ; 



(59) 



where Tfc e M(2ifc-i)x(2Lfc-i) ^j^^j RaL^-i gj^gj^ 
by the matrix partition in ( l60t shown at the bottom of the 
page. Substituting (ISSl l and (|59] | into (|57] ). we obtain 



Xk 



47r2/3- 



(61) 



which is independent of all the amplitudes. 

Note that < Xfc !i 1- Xfe is nonnegative since it is a 
quadratic form and is a positive semi-definite FIM (hence 
is T j^^); and Xk < 1 since the contribution from each anchor 
to the EFIM in {55[ is nonnegative. □ 



B. Proof of Corollary [7] 

Proof: This scenario can be thought of as a special case 
of Theorem [T] with Lk = 1, i.e., the first contiguous-cluster 
contains only one path. In this case, ( |6TT i becomes 



Xk 



1 i?2(l,l) 



47r2/32 i?,(l,l) 



Since waveform s{t) is continuous and time-limited in realistic 
cases, we have 



4(1,1) = |-i?,(T) 



= 0. 



T=0 



implying that Xk = 0, which leads to ( fT9] l. 



□ 



i?s(i,i) ti 



tk 



Tk 



i?s(l,l) 

4(1,1) 

4(1,2) 
4(1,2) 

4(1,4) 
4(1,4) 



i?.(i,i) 

i?.(l,l) 
-4(1,2) 

-4(1,4) 

i?.(l,4) 



4(1,2) 
-4(1,2) 



4(1,2) 
i?s(l,2) 



4(1,4) 
-4(1,4) 



4(1,4) 

Rs{l,Lk) 



Rs{Lk,Lk) Rs{Lk,Lk) 
Rs{Lk,Lk) Rs{Lk,Lk) 
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C. Proof of Theorem \2\ 

Proof: When a priori channel knowledge of the channel 
is available, the FIM is 

TnlAnlT^l + TlAlTJ TnlAnl 
AnlT^l A^ 



1 



k-NL 



where Anl = Ee {Anl} = diag *2^- ■ • , *m} and 
Al = Eg {Al} = diag {*M+i, *M+2, . . . , *ArJ. The FIM 
Jg can be partitioned as ( fT2] i. where A is given by ( |62] | shown 
at the bottom of the page, and 



and 



G7\/+i\I'7\/+i- 




„2wA/+l 







7% 



diag I * 



Af+l 



^2 — Af+1 



A'b 



2w^b 1 



Apply the notion of EFI, and we have the 2x2 EFIM, after 
some algebra, given by ( |63T l at the bottom of the page. From 
(|9]l, we can rewrite Hp p and Hp ^ in (fTTT l using chain rule as 



and H^ = qfe H^^ , 



(64) 



where . 



and Hrf 



Fe{Kk\dk;dk,dk) 
Fg{Kk\dk ',dk, Kk). Substituting ( l64l i into ( I63t leads to 
(ISTT i. where Xk is given by (165 al l and ( I65bb for LOS signals 
and NLOS signals, respectively, shown at the bottom of the 
page. □ 

D. Proof of Corollary \2\ 

Proof: We first show that the a priori channel knowledge 
increases the RII. Consider Xk in ( I65ab . Let 



1 



iJ^fclfc + c^S*^ 



fcifc ■ 



d,d 
~kT 



Dfc*fcDi + c^H 



,27rfc 



and 



1 



We have F^: )^ E^, since 

Ffc - Efe = 



- fcT 



p.K 



Fe(/«feMfe;0fc,6»fe) h 0, 



where 0^ = 'y dk \ ■ Hence we have Xk = 
l/[Ffci]i,i > l/[Ej:i]i,i, where [FI^Va equals (O. This 
implies that the a priori channel knowledge can increase the 
RII. 

We next show that the RIIs in ( I65al ) and (I65bl l reduce to 
( fT6b and zero, respectively, in the absence of a priori channel 
knowledge. 

When a priori channel knowledge is unavailable, H^ 
Hp ^, and H p p al l equal zero, and the corresponding RII Xk 
in (I65ab and (I65bb becomes 



*fcDT) (Dfe*fcDT) '(Dfc*felfc)} 



87r2/32SNRW 



^(1-X.)SNR«, 



for fc e A/l, and 



for k eJ\fw 



1 



□ 



A^ 



„2c.fc 



T -X- r2wl 



P,P 

^T 



Gi*iD{ 



T I „2rrl 



Gm^mD^.+c^H*^ 



Gm*mDT,+c2H^, 



Dm*mD];, + c2h;^^, 



(62) 



[ /ceA4 fceA^L 



+ c 



„2 — fe 

■'PjK 



feeA/Ni 



(63) 



l{ll*felfe + c^S^ , - (iJ*,.DT + c^H^ J (Dfc*,Dl + c^H^^ J "'(iI*.DI + c^H^^)^ }, fc £ A/L (65a) 



Afe — 



^{ll^kh + c^^id - (}l^k + c^Hl.) + c^H^,,) "'(ij^fc + c^H^,,)^ }, fc G A-NL (65b) 
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E. Proof of Corollary |5] 

Proof: The block matrices ^ and ^ in ( ITTT i for 
NLOS signals can be written as 



Vfc a« 



and H* 



where Vfc, „ G M^i^-i^ and H^_^ G TR(2i"=-i)x(2i'=-i) 
Note that corresponds to the Fisher information of 6^, 
When the a priori knowledge of h^j^^ 
gbibk^) (^{bk^)^ we claim that 

-1 



(1) 
fe ■ 

oes to oo. I.e., 



lim 



k + c 



,2 — 








0^ 



(66) 



To show this, we partition as 

,■2 



and then the left-hand-side of (l66Tl becomes 



LHS = lim 



where 



B ^ 
and 



= lim 



u'l + c2i2 + C^vl 

kfc + c^vfe + c^H^^, 

A B 
BT C 



H -1 



(kfe +c2vfc)^ (*fc 



■t2 



fe + C Vfc 



.2^.« 



1 -1 



MI + C't' 



(kfc + c^Vfc) (kfc 



When fel^'' is known, i.e., t"' oo, we have \imt2_^^ A 



0, limj: 



, B = 0, and limj 



Notice that 'S'fe = Dfe'S'feD]!. Hence, we proved our claim in 
Substituting ( l66b into ( I65bl i. we have 
lim Afc = ^(ll*felfe + c^S^ - (ll^kBl 

X (Dfc*fcDT + c2<^)"'(lT*,DT 



for fc G A/nl, which agrees with the RII of LOS signals 
in Hence LOS signals are equivalent to NLOS with 

infinite a priori knowledge of 6^,^^ for localization. 



□ 



F. Proof of Proposition \2\ 

Proof: Note that qfe, p, H^^, and H^^ are 

functions of p when a priori knowledge of the agent's position 
is available. Hence we need to take expectation of them over p 
in ([Tol l. After some algebra, we have the EFIM for the agent's 
position as ( |67] | shown at the bottom of the page. 

When the condition in ( |22] | is satisfied for the 
functions .g(p)'s: 1) qfeS^^q];, 2) qfc l]! U q]!, 3) 

qk (ll*fe + c2 and 4) *fc + H^ we can approximate 

the expectation of each function over p in ( |67| i by the function 
value at the expected position p. Hence the EFIM in ( |67] | can 
be expressed as 



Je(p) = 



fe6A4 I 



cik[Vk^k + c's',,J (*fe + H 



I ^2 ^fe \ „T 
k + C ^d,d) Ik 



— fe 



fe + c'sL) qj 



= 3p + ^ Afe Jr(<?i'fc) , 

where 0fe is the angle from anchor k to p, and Afe is given by 
(|68] | shown at the bottom of the page. Note that all functions 
are evaluated at p. □ 



-'Note that the size of ^ and SJJ ^ for LOS signals and NLOS signals 
are different for the same L^. Indeed, 3^ ^ and ^ are not associated 

with b^j^\ and hence they are in the same form as Sj^ ^ and 3^ ^ for LOS 
signals in j65at . 



Je(p)=Hp+ ^ |Ep{q,S^,,qT} + lEp{qfelT*,lfeqT} 

- 1 Ep {q, (ll*, + c2 H;;^^) } Ep + H^:,,}"'Ep | (ij* 



fe + c'HL) ql 



(67) 



Afe — 



1, -I- r2 "'^ - 



\l^k + c^Ht) (*fe + H^..) "'(ll *fe + Hi.)' } 



(68) 
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Appendix D 
Proofs of the Results in SectionHvI 

A. Proof of Theorem \3\ 

Note that this proof also incorporates the a priori channel 
knowledge. In the absence of this knowledge, the correspond- 
ing results can be obtained by removing Jp that characterizes 
the a priori channel knowledge. 

Since p and ip are deterministic but unknown, the joint 
likelihood function of the random vectors r and 9 can be 
written as 

f{v,e) = f{v\e)f{e)= [] f{^nk\e)fM\y,^). 

neAC feeA^b 

Note that /(/«„fc|p, (^) = f{nnk\dnk), and the FIM Jp from 
f{6) can be expressed as ( |69] l shown at the bottom of the 

p,<p ~ Hnk ^d,d ""nk. 



page^ where H^J^p = q„fcS3^^,qT„ H'"^- 
and = hl^ in which 



F6)(rnfc|^ ; dnki dnk) ■ 



Block matrices Hp,n, „, and H„ correspond to the nth 
antenna in the array, and they can be further decomposed into 
block matrices corresponding to each anchor: 



and 



r — n4 



Where H;;*^, 



Kk'Sd%, and 

nk 



= Fe(K„fe|p, (^;K„fc,K„fe), in which H 

Similar to the proof of Theorem |2] in Appendix IC-CI the 
FIM from observation can be obtained as dTOb shown at the 
bottom of the page, where 



G„ — [ q„4l^ ;^ 'in,2^n^2 



and 



diag{*^ 



I ^ r 



I ^ r 



} 



correspond to the nth antenna as defined in (|44] |. 

The overall FIM Je is the sum of (|69] | and ( TTOI i. By applying 
the notion of EFI, we have the 3x3 EFIM for the position 
and the orientation as follows 



E E 



\ T 
An/c Qnfc q^A; 

'^nfc ^nfc Qnfc 



^nk ^nk Qnfc 



nk 



(71) 



where A„fc is given by (|72] | shown at the bottom of the page. 

Note that in the absence of a priori channel knowledge, 
the above result is still valid, with the RJI of (|72] | degen- 
erating to ( l73T l shown at the bottom of the page, where 

D„fc = [ l2L„,-l ]. 

B. Proof of Proposition |5] 

Since qq^ is always positive semi-definite, we need to 
simply prove that there exists a unique p* such that q* = 0. 



'keK 



^nk 

"^P:P 

^nk ^ 



■ikeAfb "p,i/5 
— T 



— T 



i=p,l 

3i 



(69) 



— ..GA:, '^"^^"W I^neA/:, ^n^n'^n GiAi 

E„GA^, h„A,„GT J2neK h«A„hT hiAi 
AiGj Aihj Ai 



Gat^Ajv, 



(70) 



(72) 



'iL - i^nk'Dlk) (D„fc*„fcDT,) '(D„fc*„fe) } l„fc/c2 , LOS 

, NLOS 



(73) 
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Proof: Let p be an arbitrary reference point, and Therefore, we can verify that the EFI for the orientation in 



P* =P + gM: 



where g((^) = [gx{^) gy{v)Y ^ and gx{ip) and gy{ip) denote 

the relative distance in x and y directions, respectively. Then, ^2 * ~ 

hnk corresponding to p can be written as a sum of two parts ~ 2^ 2^ "'^ "^^ ~^ 2^ 2^ "'^ "'^ 

hnk^Kk+hnk, =J*{(p). (76) 

where /i* corresponds to p* This shows that the EFI for the orientation is independent of 

^ ^ the reference point, and thus is the SOEB. 

/i*^ = — Ax„(p*, 1^) coscpnk + —Ayn{p*,(p) sincpnk , We next derive the SPEB for any reference point given in 



dip dip 



The 3x3 EFIM in (|7B can be written, using (|74| and 



and (|76l), as 

hnk = '^9xiv) COS(t>nk + '^9v{v) sin Je(p, (/?) 

;.T 



J q 



- coscpnk + ffy Sill 0„fc - g q„fc . Using the equation of Shur's complement 1661, we have 

Hence, q corresponding to the reference position p is given 1 /,i^x/^i>xT 

by J"(P)=J 

q= Z ^^nkKk'^nk+^ ^Kk hnk Hnk, (74) =J"^ + — ^gg'^. (77) 



nGA4 k£M, n£Mi k£jVb 



Since the translation g((p) can be represented as 



= q =q 
and q can be written as g('p) = l|P ^ P* 



cos(v5 + ipo) 
sin((^ + ipo) 



^ Z/ Z/ *lnfc g -^"fc q^fc where t/So is a constant angle, we have ||g| 



P P 



Then, by taking the trace of both sides of (ITTI i. we obtain 



• T • 

neM,fceA4 neM, ^(p) = 7'(p*) ' ^ ^ 



Je((/5) 

Since 2]„eM, Je,n ^ 0, we have q* = if and only if ^ ^^p*^ ^ Hp _ p*||2 . ^^^^ 

g = Z "^<='" ) ^' 



implying that there exists only one g, and hence only one 

g(^), such that q* = 0. Therefore, the orientation center p* Take the array center po as the reference point. 



□ 

D. Proof of Proposition \5\ 

Proof: T 
and we have 



IS unique. □ 

Z ^"fc ^ Z 3~^^n(P0'V') COS(/)„fc 

C. Proof of Corollary\5\ «eM, neM, ^ 

Proof: We first prove that the SOEB is independent of + ^ — A?/„(po, (ys) sin 4>nk 

the reference point p. It is equivalent to show that the EFI for neA/; 
the orientation given by ( |27] i equals the EFI for the orientation d / \ 

based on p*, given by ^ I Z ^^"(PO'*!^) ) cos0„fe 



neM,fceA^b + A I ^ A?;„(po,(^)) sine 



'rifc 



Let J = E neA/;Je,n- From d?!]) and (|75]l, we have q = 
q = J g, and hence ~ ^ ■ 



q^J ^q = q^J ^q = q^ g = E E ' 



Consequently, 



On the other hand, we also have 



q = E E ^khnk qfe = E E '^"'^ j Afc qfe = , 

Z] E/ ^nkhnkhnk q g = . implying po — p*, i.e., the array center is the orientation 

n6A/;, fceTVb center □ 
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Appendix E 
Proofs of the Results in Section[V] 

A. Proof of Theorem |4] 

In the presence of a time offset, the FIM can be written as 
dTSl l shown at the bottom of the page, where 






















■ QT 


Jp = 


— 1 T 
' — ^p,*^ 











W^b T 










Applying the notion of EFI, we obtain the 3x3 EFIM 



Je(p,S) 



EfceAAb qfc ql EfceAAb m 



keAfi ^k 



Afc + Se 



where Xk is given by ( I65bl i. and another step of EFI leads to 
(m and (O. 



B. Proof of Theorem \5\ 

We consider orientation-unaware case, whereas orientation- 
aware case is a special case with a reduced parameter set. The 
FIM using an antenna array can be written as ( |79] | shown at the 
bottom of the page, where 1„ 1^ ^ 1^2 ••• ^1,nA^ 
and Jp is given by dSOl l shown at the bottom of the page. 
Applying the notion of EFI to Je, we obtain the 4x4 EFIM 
in ((371). 



C. Proof of Corollary |S] 

We incorporate the a priori knowledge of the array 
center and orientation into ( |37] |. and obtain the EFIM 
in far-field scenarios as dSTl i shown at the bottom of 
the page. Recall that in far-field scenarios, po = 
p*, implying that Y.n<^N■,Y.ke^^b^nkKk^.nk = and 
LneA/:, LfeeA^ ^nkKk = 0. Also, we have A„fe and 
0n/c = 0fe for all n, and hence the EFIM can be written as 
( l82b shown at the top of the next page, where hnk and is 
a function of po. 



iT 



Gfc*fclfe Gi*i 

T,T,.l. IT^^ 
*lll #1 



l]vb*^b 



+ Jn 



(78) 



JVb 









SnGA/; G„A„1^ 


GiAi • 


GjVaAjv^ 




SnGA4 h„A„G!^ 


SnGA4 h„A„h^ 


X^nGA/; hnA„l^ 


hiAi • 




1 




X^nGA^, lnA„h^ 


X^riGA/:, InA-nl^^ 


liAi • 


IwaAjv,, 




AiGl 


AihT 


Aill 


Ai 














A-w., 



(79) 



Jp = 



SriGA/;, SfcGA/i, 



feGAb "P,P 
ynfe T 
"■p.V 



— T 

"pa 



— T 

"p,Mi 



/cGA4 



nGA^, Z^fcGA/b 



— T 



— T 











■^pa 



"P,Wa 
■^V.Wa 



^A^a 



(80) 



X^TiGA^a SfcGTVb ^nk^nk^ 'nk + "P SriGA^a SfeGA/i '^nkhnklnk EiiGA^a SfeGA/i -^nfeq^fe 



jAnay-B 



SnGTVa X^feGA/i ^rikhnk^!nk 
2nGA4 SfcGA^b ^"fc^n* 



, X^feGA/i ^nkhnk 



J2neK YkeAfb '^■'ik^nk + "V J2neM., 

SnGA4 SfcGA/b ^nkhnk 2nGA4 SfeGA/L ''^"^ + 



(81) 
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jAiray B ^ qT Y^neK YkeNi, ^kh'lk +■='¥' ^ 



(82) 
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